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Abstract. In this paper we investigate the asymptotic behavior 
of the semi-classical limit of Wigner measures defined on the tan- 
^ , gent bundle of the one-dimensional torus. In particular we show 

^**^ | the convergence of Wigner measures to the Mather measure on 

the tangent bundle, for energy levels above the minimum of the 
effective Hamiltonian. 

The Wigner measures y,h we consider are associated to ip^ , a dis- 
tinguished critical solution of the Evans' quantum action given by 
^ = ah e'*, with a h (x) = e <M ^\ u h ( X ) = p.*+<M+^M 
and Vh , satisfying the equations 

oS: -^+ 1 -\P + Dv h \ 2 + V = H h (P), 

en; f ^ + ^ P + Dvl f + v = H h (p), 

where the constant Hh{P) is the h effective potential and x is on 
the torus. L. C. Evans considered limit measures \iph\ 2 in T™, when 
h — > 0, for any n > 1. 

We consider the limit measures on the phase space T n x R™, for 
, n = 1, and, in addition, we obtain rigorous asymptotic expansions 

?-h ' for the functions Vh, and v^, when h — > 0. 
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1. Introduction 

Consider the mechanical system defined by a Lagrangian of the form 
L(x, v) = \\v\ 2 —V(x)—Pv, with x in the n-dimensional torus TP = 
v G M. n , and fixed P G M n , and the associated Hamiltonian given by 
H(p,x) = ~(P + p) 2 + V(x). Here, V stands for the potential energy, 
which we assume to be symmetric in order to simplify our arguments. 

Before we present our results, let us review some basic facts of the 
Aubry-Mather theory. Let Ai denote the set of probability measures 
on the Borel a-algebra of T n x W 1 . 

The Mather problem on T n (see |Mat] . |CI] and |Faj ) consists in de- 
termining the probability measures fi G M. which minimize the action 



:i) / 

</T n x] 



1 \v\ 2 -V(x) 



d/i(x, v) 



2' 

among the probabilities /ieM such that (the holonomic condition) 



J v.D(f>dfi = Q, V0 G C^TT) 



and, 

v d\i = V, 



for a fixed vector V. 

For dimension n > 1, these minimizing measures on the tangent 
bundle, called Mather measures, do not need to be unique and are 
supported on sets which are not attractors for the flow. They are 
invariant by the correspondent Euler-Lagrange flow. 

It is known that the Mather measures on T n x M n are supported on 
graphs (the projection on T n is injective). The probabilities we get 
when we project on T n are called projected Mather measures. They 
are probabilities on T n . In most of the cases the projected probability 
is not absolutely continuous with respect to the Lebesgue measure on 
the torus. 

The Mather measure described above can be alternatively obtained 
as a minimization of the action of E-L invariant measures without con- 
strains considering a new Lagrangian of the form (see |CIj ) 

L(x,v) = |H 2 -V{x)-Pv. 
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The result presented by L. C. Evans in jEvlj considers the projected 
Mather measure, but his methods do not extend to the analysis which 
we will consider here on the tangent bundle. 

An example of Mather measure on the tangent bundle can be shown 
on figure 1. In this case consider the only E-L invariant measure with 
support on the top curve in figure 1. The periodic curve which is 
inside the open set determined by the two separatrices does not carry 
a Mather measure. 

In the above case the projected Mather measure can be computed 
by action angle variables (see section 7.3 page 110 (7.24) in |PRj . or, 
section 7.2 page 159 (7.22) in jOA] l 

Now consider the Hamiltonian H associated to L, that is given by 

H (p, x) = sup[pv — L(x, v)], 

and its associated Hamilton- Jacobi equation 

H(V(f)(x),x) = c, 

where c is a constant. In general, we cannot find C 1 solutions of this 
equation and, therefore, we consider viscosity solutions, as follows. 

We say that a function : U — > R is a viscosity sub-solution in U 
of the equation H(V<j)(x),x) = c if, given any C 1 function rj : U — > R, 
and any x G U, we have 

H(Vrj(x ),x ) < c, 

whenever <f> — r\ attains a maximum value at Xq. Dually, we say that 
a function : U — > R is a viscosity super-solution in U of the same 
Hamilton- Jacobi equation if, given any C l function rj : U — > R, and any 
xq G U, we have H (Vr](xo), xq) > c, whenever — t] attains a minimum 
value at xq. Finally, a function is a backward viscosity solution of the 
Hamilton- Jacobi equation if it is both a viscosity sub and super-solution 
(See pM], [Ea], PS2], [EG], and [A~T]). 

Now let us define a forward viscosity solution of i7(V0(x),x) = c, 
see [Fa] . To do that we need to define 

H(p,x) = sup[— pv — L(x, v)]. 

V 

Let be a backward viscosity solution of H(V<f)(x),x) = c, then we 
say that 0* = —0 is a forward viscosity solution of H(V(p(x),x) = c. 
Note that H(p,x) = H(—p,x). 
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L.C. Evans has shown that for a Lagrangian of the form L(x,v) = 

^ V{x) — Pv, the solutions Vh and v* h of (J2J) and (j3J) converge (maybe 

in subsequence), respectively, to backward and forward viscosity solu- 
tions, (ft and </>*, of H(V<p(x),x) = H(P), where H(P) is sometimes 
called the Mane critical value, see [Evl] . 

If (j) : U — > R is a viscosity solution of H(V(f>(x),x) = H(P) then 
the equality H(V<f>(x),x) = H(P) holds, in the classical sense, at every 
point x of differentiability of 0. (See Corollary 4.4.13 in |Faj.) 

In the case the Mather measure is unique, the viscosity solution is 
unique, up to an additive constant (See Sections 7 and 8 in [Fa], Section 
4.9 [CI], [Evl] . [BG] . and |Gom3j ). and the solution v h of © does 
converge to the backward viscosity solution of H(V<f)(x),x) = H(P), 
as h — > 0. 

There exist a value P cr u such that for P > P cr it there exists viscosity 
solutions (see [Fa] ) on the corresponding level of energy H(P). For the 
case of a potential with just a maximum in the one- dimensional torus 
the constant energy lines are described in figure 1. The periodic curve 
in the top corresponds to P above the critical value P cr u- The periodic 
curve which is inside the open set determined by the two separatrices 
not. 

We will consider here only the one-dimensional (that is, n — 1) case 
for a certain general class of potentials for which the Mather measure 
is unique, and the projected Mather measure is absolutely continuous 
with respect to the invariant measure. In this one-dimensional setting 
the Mather measure is the unique measure on a certain energy level 
H(P) (depending on P) which is invariant by the Euler-Lagrange flow. 

For level of energy H(P) which projects on the all one-dimensional 
torus one can consider the function p + (x) which graph determines the 
curve on the tangent bundle. The integral §_y 2 p + (x)dx = P, relates 
the energy level H(P) to P. The function H(P) is differentiate on P 
[EGT] . 

We point out that our result does cover a large class of interesting 
cases, such as the one-dimensional Lagrangians L(x, v) = \ \v | 2 —V(x) + 
w x (v), with w a nontrivial closed form, for which the projected Mather 
measure is absolutely continuous. 

Evans introduced a quantum minimization problem for his quan- 
tum action which is similar to the classical minimization problem of 
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Mather's theory. He considered the semi-classical problem of limit 
measures on the configuration space TP, showing that in the semi- 
classical limit with h — > 0, when the probability associated to the wave 
iph converges, then the limit measure is a projected Mather measure 
[EvTj . 

To describe our main result we need to recall some facts from |Evlj . 
We analise the Wigner measure associated to the critical function 

. Px+v{x) 

i) h {x) = a[x) e l h 
of the Evans's quantum action 

Sty] = [ \h 2 \Da\ 2 +[\\P + D x v\ 2 - V{x)] a 2 dx, 

under the constraints 

/ a 2 dx = 1, 

stationary current 

div(a 2 (P + D x v)) = 0, 
and total current intensity 



/ a 2 (P + D x v)dx = Q. 



The function a above take only real values. 

The analysis of the properties of the critical solutions described above 
is the so called Evan's quantum action minimization problem. 

The Evans' critical solution (see |Evl] ) is given by the periodic func- 
tion i/j h = a,he l ~h, where dh(x) = e~ 2?r~ , and Uh(x) = P ■ x + 
v h {x)+v h {x) ^ g ere ^ VhjV * satisfy the equations 

(2) - ^+ l -\P + Dv h \ 2 + V = H h (P), 

(3) ^ + l\P + Dv* h \ 2 + V = H h (P), 

where the constant Hh(P) is called the effective Hamiltonian. 
We also assume that — Vh is normalized, that is, 

r v* h ( X )-v h ( X ) 

(4) e h dx = 1, 



T" 



and, therefore 



1^(^)1 dx = I. 
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In the present case we consider the limit function vq which is a vis- 
cosity solution (also almost everywhere differentiable) and satisfies div 
((P + v'q) do) = 0, where a is the projected Mather measure (see |Evl] 
section 3, |BG] or |EG1] ). We will comment more about this fact in 
the last section. 

The function ip h is not periodic. 

N. Anantharaman |Alj . |A3] and L.C Evans |Evl] had previously 
(and in arbitrary dimension) studied the convergence on the configu- 
ration space of the probability measures \iph(x)\ 2 (x) dx. 

One of the main motivations in the study of the critical solution given 

u / ;i .,, , x , . v*(x)+v h (x) 

by iph = dh e h , with (Xh{x) = e ^ } Uh[x) = P ■ x + — 2 — — , as 
h — y 0, is its relation to the viscosity solutions of the Hamilton- Jacobi 
equation which was described before. 

C. Evans [Evlj considered probabilities "on the n-dimensional torus". 
We are interested here in the semi-classical limit of the Wigner mea- 
sure (to be defined soon) associated to iph, when h — > 0, and its re- 
lation to the Aubry-Mather measure "on the tangent bundle" of the 
one-dimensional torus. (See ( |Mat j . [CI] . [Fa) , and |BGj .) This means 
that we consider the problem on the phase space T n x R n . 

We point out that the quantum minimal action problem is closely 
related to the Bohm's viewpoint of Quantum Mechanics, which is con- 
nected to the classical Hamilton- Jacobi equation with the addition of 
an extra diffusion term. (See Section 19 in [Jo], |Re] and |Voj .) Observe 
also that we are considering an action which is quite different from the 
action used by F. Guerra and L. Morato in [GMJ. 

Before we proceed, we discuss some background material on the semi- 
classical analysis and Wigner measures on the torus T n = M n /Z n , which 
we identify with [— |, |)™. 

Given a state ip G L 2 (T n ), one would like to compute the averages of 
observables, such as momentum or energy. Here, observables are linear 
operators b and their averages are given by {ip, Injj). 

Frequently, quantum observables are pseudo-differential operators 
associated to a smooth symbol a(x,p). In the periodic setting, the 
momenta p are quantized and, instead of taking values in IR n , they 
belong to the lattice KL n + In the pseudo-differential calculus, there 
is not a unique way to associate an operator to a given symbol. A 
common choice is to associate to a symbol a the Weyl operator a w 
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given by the Weyl quantization rule 

P «/T' 1 ^T" 



which has the advantage of associating self-adjoint operators to real 
symbols. (See |dGo] . and |A3] for related results.) Notice the mul- 
tiplicative factor 2 ir of the momentum variable, due to the product 
coordinates of the torus. In our one-dimensional setting, 



i i 



P £hz+£ * 



dx. 



If a is simply a function of x, then 

(ip,a w t/j) = / a(x)\ip\ 2 dx. 



and if a depends only on p, the average value is 
where 



pehZ n + P_ 



27rip-x 

V'(p) = / ip{x)e ~ dx 



are the Fourier coefficients. 
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Given any ip such that j \ijj\ 2 dx = 1, we define on T n x (h Z n + ^) 

by 



Note that < ip, a ip > gives real values when applied to smooth real 
functions with compact support (See |Zhj ). The value corresponding 
to a nonnegative real function a it is not necessarily nonnegative. 
It is well known that if f Jn \ip\ 2 dx = 1 for some if) G L 2 (T n ), then 



For fixed ip it is known that limit points, when h 0, of Wigner 
distributions are positive measures |Zhj . 

We denote by fi[ h the Wigner distribution associated to ijjh on T n x 
MJ 1 given by 



We want to consider for each value h the Wigner distribution Hh = 
fij^ h associated to the Evans' critical iph- Note that the normalization 
we considered before for — does not create indeterminacy for W^ h . 

A very complete description of what is known when h — > 0, in prob- 
lems related to the projected Aubry-Mather measures is [Alj . |A3j . 

In |AIPS] the case where do exist several different Mather measures 
is analyzed and the question of selection of the limit projected Mather 
measure is considered. 

In |Evlj the problem is described in a different setting and in terms 
of the quantum action. This result is for the torus of dimension n. 

In the general case the Mather measures are singular with respect to 
Lebesgue measure. 




It follows that 
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On the geodesic case for negative curvature manifolds these measures 
are supported on geodesic laminations (appendix |A2j ) and a Large De- 
viation Principle is known under the hypothesis of uniqueness of the 
maximizing probability. For nonnegative curvature the Large Devia- 
tion problem around a closed geodesic with only points of curvature 
zero is analyzed in |LRj . 

Our results are all in the one- dimensional setting, and our main result 
is the following. 

Theorem 1. Let L:TxI->lko Lagrangian of the form 

L(x,v) = \\v\ 2 -V(x)-Pv, 

with fixed P > P cr u, and a smooth potential energy V symmetric with 
respect to the origin that has a unique non- degenerate minimum point 
at the origin. 

If fih denotes the Wigner distribution associated to the Evans ' critical 
iph, then fih converges to the Mather measure on the tangent bundle. In 
this case, the Mather measure is the unique invariant measure for the 
Euler-Lagrange flow on a certain energy level (which depends on P). 

Moreover, consider a smooth function f supported on a compact set 
AcTxi 

(1) If the support of f does intersect the energy level H(P), then 



2 

f(x,p) dfi h (x,p) 2n^2 J 



f(x i (p) J 2irp)p + (x) d 

\(p+)'(Xi(p))\ P + (Xi(p)) 



for certain functions p + and Xi(p) with p + {xi(p)) = 2irp. 

The value x appears in a natural way and is responsible for 
the normalization which assures us that we get a probability in 
the limit. The right hand side is just the canonical expression of 
the invariant density for the Hamiltonian flow via action-angle 
variables. 

(2) If the support of f does not intersect the energy level H(P), 
then 

f(x,p) d/i h (x,p) 0. 

In both cases case we are able to estimate the speed of convergence. 
Our measure convergence claim is in the following sense: fih con- 
verges to /i if J fdjih —> J fdfj,, for each smooth functions f with the 
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property that its support does not intersect two certain lines p = p^, 
and p = Pmax t° be defined in the following section. 



The symmetry hypothesis on V simplifies our computations, but our 
method can be adapted to more general situations (taking up much 
more pages). For instance, just below (1151) this symmetry hypothesis 
results in the same property for the viscosity solution, and this makes 
the analysis of the problem much more simple. 

In the proof of our theorem we obtain an estimate of the velocity of 
the convergence of J f(x,p) dfih(x,p) for a certain class of test function 
f(x,p) with compact support, where /i^ is the Wigner distribution as- 
sociated to the family ip^. In addition, we also establish the asymptotic 
behavior of the solutions of ([2]) and ([3]) when h — > 0. This requires a 
rigorous asymptotic expansion analysis which, to the best of our knowl- 
edge, cannot be found in the existing literature. 

It is true that in this one-dimensional case we could avoid the men- 
tion of Mather measures, and consider only invariant measures on con- 
stant energy levels but since we are going to use the setting and the 
results of |Evlj . which are presented in the context of the Aubry-Mather 
theory, we believe that our discussion so far could be helpful. 

Note that if one knew in advance that the limit of the Wigner dis- 
tribution is an invariant measure for the hamiltonian flow, the result 
would be trivial, because we are in a situation where there is a unique 
invariant measure on each energy level. However, h is only a 0(h)- 
quasimode of the quantum hamiltonian (this was proved by C. Evans 
[Evlj ) and we do not know a priori that the limits of the Wigner mea- 
sures are invariant. 

The proof relies on the facts that the phase space picture is com- 
pletely explicit, that the Aubry-Mather measure is unique, and that 
the solution to the Hamilton- Jacobi equation is unique, smooth, and 
explicit; thus, this approach is specific to dimension 1. 

The semiclassical limit of Wigner measures that minimize the quan- 
tum action was considered in [GVj . This was inspired by the previous 
work jEvl] which only uses classical tools. We now revisit this prob- 
lem and compute the rates of convergence (above the critical level of 
energy) for one-dimensional Hamiltonians. We will consider P > P cr u, 
where P crit = inf{P : ~H(P) > minI7} > 0. 
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Concerning possible extensions of our results to higher dimensions, 
we acknowledge that the control of the stationary phase method in the 
first part of our work, although being much more difficult, could pos- 
sibly be obtained. Anyway, in the last part of our work (see Section 
3) we have to control an asymptotic expansion, for which we need a 
Mather measure that projects on the entire torus. We strongly be- 
lieve that this control is not obtainable with our methods without that 
assumption. An extension in this direction will require a completely 
different approach. 

This paper is organized as follows. In the next section we derive 
the correct asymptotic behavior of the integrals over Wigner measures 
when h — > 0. In the last section we describe the asymptotic expansion 
of the function Vh in terms of h, which we will need for the proof of 
Theorem [U 

2. The asymptotic limit of critical Wigner measures 

The critical Wigner measure is the Wigner measure associated to 
a critical solution ip^. For simplicity we will omit the word critical 
from now on. We want to investigate the asymptotic limit of Wigner 
measures, when h — > 0. 

Remember that W% h (x,p) = L^(x + — |)e f^dy, hence 

for any continuous / with support on a compact set A C T x R 

"fafo+f )~ v h<- x +% )+ v h( x ~% )~ v h ( - x -j ) 
2h 



I f(x,p)dfi h (x,p) = ^2 f( x i 27l P)e 

" r-n P " ^ " T 



■e 1 * + 2h J e h dxdy. 

Note that, as p E {hm + m E Z}, we can write 



/ f(x,p)d^ h (x,p) = ^2 f{x,2irp)e 



p&hZ+- p 



v h ( - x+ il- v h( x +i)+ v i( x -^)- v h( x -^) 



. vHx+^)+v h [x+y)-vi<.x-y)-v h (x-y) 
■e 1 \ e zmm v dxdy. 

Note that the integrand is periodic on the variable y 
We point out that the analysis of the case of integrals of functions 
depending only on x, namely, f(x,p) = f(x), is well known. It follows 
from the results in |A1] and |Evlj . In particular, the projection of the 
above Wigner measures converges to the projected Mather measure (on 
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the torus) |Evl] and large deviations are also well understood (in the 
case the Mather probability is unique) according to Proposition 3.11 

In the one-dimensional case (above the x-axis) the level of energy 
consist exactly of one periodic trajectory when P > P c nt- Indeed, note 
that if the energy is high the modulus of the velocities on this level 
of energy is high, and, therefore, by conservation of energy the level 
of energy has to be a curve like the one on the top of figure 1. So, it 
follows from the graph property |CI] [Faj [BG] [Gom3j that the Mather 
probability has support in this periodic trajectory. 

The case of P G (—Peru, Pent) (where P attains the minimum value 
of H) requires a different analysis and will not be considered here. 

We will show later that W^ h converges to the Mather measure which 
has support in such trajectory. 

For P > Pent, let us find the backward viscosity solution of H(<fr'(x),x) 



H(P). If we denote by p+(x) = J2(H(P) - V{x)) > then 

f x 1 

4>p{x) = (j)(x) = / ^ pp(s)ds - P (x+ -) 

is a solution, and hence a viscosity solution, of H((j)'(x), x) = H(P). It 
is easy to see that in this case the forward viscosity solution <fi* is such 
that 4>* = cf>. 

The function is such that its graph determines a level of energy 

H(P)(a,s for example the curve shown on the top of figure 1). 

i 

Note that Pp{x) = P + <f)' p (x). Then, we have j\ Pp{s) ds = P. 

Let us denote by p^ in and p^ax the numbers such that p^in — 
Pp(x) < p^ ax for all x e T, note that p^ in > 0. 

The Mather measure with rotation number Q = DpH(P) is defined 
by the unique periodic trajectory supported in the graph (x,pp(x)). 

Definition 1. Let f and g be functions of the variable h > 0. We say 
that f(h) = 0(g(h)) as h — if there exist a constant K > and 
5 > such that \f(h)\ < K\g{h)\ whenever < h < 5. 



Definition 2. Let f(h) and i/j(h) be two functions satisfying 
f 

lim — — = 1, then we say that f is asymptotic to vp, or ip is an 
h^o ip{h) 

asymptotic approximation to f , and, we write f{h) ~ ip(h), as h — > 0. 
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In the last section we will study the asymptotic behavior of the 
functions Vh and v^, the solutions of (EJ) and (EJ). 

We prove that for any integer 7 > 2 there are functions functions 
Vi,v*, i = 0, 1, and w 1 , w* such that the following functions 

(6) i'l := v + h vi + h 2 w^ , vj* :— Vq + h v{ + h 2 w*, 

approximate uniformly the functions Vh, v^, to a order depending on 7. 
More precisely, we can choose u> 7 and given by a finite sum of the 
form 

7 7 
w 7 = ^^S' w * = 

i=2 J=2 

where all functions Vj,v* are smooth. Then if we define =: Vh — 
and ^ =: — vj*, by proposition [5] of the last section, we have that 

(7) \g h {x)-g h {x h )\=0{hi), 
and 

(8) \g* h {x)-gl{x h )\=0{hl). 

Note that the error term is only 0(h?) in the uniform norm, and we 
do not make any claim on boundedness of its derivatives. 
In addition, we will also show that 

(9) 4>p = 4> = v = Vq , vi = -u*. 

Remark. By the normalization hypothesis (jlj) for Vh and v^, for each 

v h( x h^ ~ v h ( x h^ 

h, there exists Xu £ T such that e s = 1, then vl(xh) — Vh(xh) = 
0. This implies that 



(!0) g* h (xh) - 9h{xh) = vl(x h ) - vl*(x h ). 

Now we will use the notion of asymptotic approximation to simplify 
the expression of the Wigner measure. 
We define 

r r v * h(x+ y ) - Vh(x+ v )+v i (x -y ) - Vh(x -y ) 
F Pt p(h) = F(h) := f(x,2irp)e 
Jt Jt 



■e 1 ^ 2h hi dxdy 

and 

F P3 {h) = F{h) := / / /(x^TT^e^^+^^-a^^eS^^dxdy, 
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where Sp(x,y) := P- y + v (x + |) — v (x — |) — 2irpy, and x := limx h . 

We point out that the value x appears in a natural way from the 
normalization hypothesis described by Xh (see Remark above) and as- 
sures us that we get a probability in the limit. It has a multiplicative 
effect on the limit measure and therefore its influence does not depend 
on the subsequence we choose (see expression (1) in Theorem 1). 

In the following proposition we use the results obtain in the last 
section, that is, the asymptotic expansion of Vh and v^. 

Proposition 1. Fix 7 > 3. Then there exists a smooth function 77 

F(h) =F{h) + h [ [ r]{x, 27rp, h) ei-Mx+^-Mx-^+iM^^s^y)^ 



Ji j j 
+ (/^- 2 )/ 2 ). 

Furthermore, the function 77 is uniformly bounded for all h, together 
with all its derivatives. 

Proof. First note that 

v]S? + 1 ) ~ Vh(x + |) + vl(x - I ) - v h (x - I) 

2h 

= 9* h ( x + |) ~ 9h(x + |) + 9* h (x - |) - g h (x - f ) - 2g* h (x h ) + 2g h (x h ) 

2h 

*(x + I) — v7(x + %)+i>2 



Wi* + I) ~ + I ) + vl*(x - I ) - vl(x - 1 ) + 2^(x fc ) - 2g h (x h ) 



2h 

Using equations (jUJ) and (|T0|) . we get 

^(X + |) - <X ^(g - f ) 

+ |) - + |) + - |) - 0h(z - I) - 2g* h {x h ) + 2g h {x h ) 



2h 



v d x + |) - ^(^ - |) + 2vi(x h ) 



+ t|K(z + |) + w^{x - |) - w 7 (x + |) - w; 7 (a; - |) - 2w*(x ft ) + 2w 7 (x fc )]. 
Let 



w*(x + |) + w*(x — I) — U7 7 (x + f ) - Wj(x — |) — 2w*(a;/ l ) + 2w 1 (xh) 
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Using equations (J7|) and (jHJ) we get 
(11) e ^ft 

7-2 

= e ~ Vl ( X+ 2 S >~ Vl (- X ~2^ +2vi ( X h) e hL0 i( h ' x 'y) e °( h ^) 
_ e -v 1 (x+%)-v 1 (x-lt)+2v 1 (x h ) 

.y.\ — in It — y.^ 



_|_ e -vi{x+%)-vi{x-&)+2vi(x h ) e hui(h,x,y) _ -Q 

7-2 

_|_ e -«l(x+f )-«i(a;-|)+2i)i(x h ) e /iwi(/i,x,3/)^ e O(/i^ _ ) _ 



7 — 2 

We should note that the last term of the previous equality is 0(h^~) 
in L°°. The second term is 0(h) in the C k topology for any k. 
By the same kind of reasoning we get 

v* h (x + |) + v h (x + |) - <(x - |) - - |) 



(12) 



2h 

vq(x + |) - ^o(s - |) 



2 

+ 

+ 



+ |) - W *( X ~ |) + ^7^ + f ) - ™7( X - I 

gjtfr + I) - gfcM - 9h(x - |) + gftK) 
2/i 

ffhfo + |) ~ 9h(x h ) ~ ff/>(x - %) + ghfo/Q 
2h 



Define 



. w*(x + $) - w*(x - *) + w^(x + V) - w^{x - *) 
uj 2 (x, y, h) = — £ ! . 

Then 

7-2 

_ e iSp{x,y) _|_ e iS$(x,y) — 1) + ^g^^OfliT") _ 

7 — 2 

As before, the last term of this identity is 0(h~^~) in L°°, the second 
term is 0(h) in the C fc topology for any k. 

Therefore, by combining (TTTj) and (TT3|) we obtain the result. 

□ 

Let us denote by 

((x,y) := gl—"! (a;+f f )+2^i (s)] _ 
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Hence 



p£hZ+ p 



I fW= E / / f{x,2irp){{x,y)et s *^dxdy. 



T J J 



We have to estimate the asymptotic of the limit \im h ^ If(h). 

In this paper we will use several asymptotic expansions for inte- 
grals, namely the non-stationary and stationary phase methods. These 
expansions are valid for smooth functions and are uniform depend on 
bounds on a finite number of derivatives of the function involved. Since 
in all cases the non- stationary or the stationary phase methods are ap- 
plied to (fixed) smooth functions all error estimates are uniform. In 
the use of oscillatory integrals it is possible to make expansions of any 
order, if the derivatives of order k of the functions involved decay fast 
enough. In our case all functions are C°°, and then one can use the 
C°° topology, or C k topology for any k large enough. 

To estimate If(h) we will use a two dimensional version of the sta- 
tionary phase method in the variables (x,y). Let us fix p e KL + 
the main contribution in the integral on x and y above is due to the 
term e^ s ^ x ' v > . We have to find the points (x, y) which are critical for 
Sp(x,y). These are the points (x,y) that are solutions of the system 



(14) ^(x, y) = P + y o (x + f ) + \v' Q {x - |) - 2 np = 0, 
and 

(15) ^(x,y)=v' (x+ y -)-v' (x- y -) = 0. 

Note that, by the definition of Vq(x) = </>(x), we have v' (x) = p + (x) — 
P. If x in the equation (IT5|) . by the symmetry of the function Vq, 
the only solution for v' (x + |) — v' (x — |) = is y = 0. And, if y ^ 
in equation ([15]) we must have x = 0. 

Now, if y = 0, the equation ( [141) becomes 

(16) p + (x) - 2np = 0. 

And, if x = 0, by the symmetry of the function p + , the equation ( [141) 
becomes 

(17) p+(|)_ 27 rp = 0. 



W-MEASURES AND THE SEMI-CLASSICAL LIMIT TO A-M MEASURE 17 



For each p let us call x±(p) and X2ip) the two points such that 
p + {xi(p)) = 2np, i = 1,2. 

Remark: Note in figure 1 (for the top curve) that for any given p 
(in the image of the projection in the y coordinate of the points in 
the corresponding level of energy) there are two points x±, x-i such that 
(x±,p), (x 2 ,p), are on this level of energy. Due to the convexity of the 
term | \p\ 2 in the Hamiltonian H(x,p) the level lines have this shape 
for the kind of symmetric potential we consider here. 

Note that p + (0) is the supremum p max of the function p + . Moreover, 

P ( o) Pmin- 

For potentials V with a larger number of points of maximum we 
could have a larger number of pre-images of p by p + but the reasoning 
would be basically the same. 

Therefore, we have the following pair of possibilities for (x, y) satis- 
fying both (ED and fTTol) : 

(1) y — and x = xi(p), or x = ^(p), which are the solutions of 

(USD, 

(2) x = and y = 2x 1 (j3), or y = 2x 2 (P), which are the solutions of 

(ED- 

Summarizing, we have two possibilities: 

(i) For p e hZ + £ such that 2np e K^pLJ, the function Sp(x, y) 
has 4 critical points: (xj(p), 0) and (0, 2xi(p)), for i = 1,2. 

(ii) For p e hZ+^ such that 2np ^ \PmtniPmax\i the function Sp(x, y) 
has no critical points. 

In order to prove the theorem[T]we need some auxiliary computations 
that we discuss now. As the limit of the Wigner distribution is a 
positive measure ( |Zhj consider quite general hypothesis) it is enough 
to analyze the following two cases: 



Case I) If the support of / is contained the strip T x (p mini Pmax] 



Case II) Suppose that the support of / is contained in T x [c, d] , with 
c, d such that [c, d] fl [p min ,p^ 



max! 



Case I) First of all we point out that when analyzing an estimate for 
the expression 

7 /W = E ii f( x > 2n p) ^ y) eisp(x,y) dxd yi 



T JT 
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we will have to consider for each h, in a separate argument, the contri- 
bution of 

/ / f{x,2irp d ) ((x,y)ei s ^dxdy, 
Jt Jt 

where pd = Pd(h) is the first point in KL + j- bellow p^ ax . We will 
estimate this term in the end of our reasoning. 
We will denote 

(18) V / / f(x,2*p)C{x,y)ei s * x >ridxdy, 

. Jt Jt 

the sum is taken over all the 2 up < pd- 

Given e > consider rj € , a C°° function such that, r] e (y) = 1 for 
\y\ < e/2, and, ^(y) = for \y\ > e. 
We write 

*i(e) = Z) / / V e (v) f(x, 2np) C(x,y)- 
. Jt Jt 

■ r p y , f(](^+i)-"0(^-f ) 27Tip V ] 

■e 1 ^ h h i dxdy, 

and 



< 



T JT 



■e 1 h ) e h dxdy. 

Lemma 1. For any fixed e, we have limine) = 0. More precisely we 
show that the expression goes like 0(h°°). 

Proof. 

Consider a C°° function g(x,p) with compact support with all its 
derivatives with respect to p uniformly bounded in h. Then, for fixed 
x and i/^0 and A a bounded set we have 

2tt i (p-\-h) y 2tt i py 

2iYipy ^ — -\ 6 h — 6 h 



e 2n i y _ I 

p£A,< P&A,< 

Eg(x,p — h) — g(x,p) 2wip V ^ 0(h) 2nij y 
g27r iy _ I Z-^/ g27r iy _ I 

p&A, < peA, < 



e h . 



We are using above the notation: J2peA < * s the sum over the P om ts 
in A which are smaller then pd. 
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Note that in the above reasoning it is fundamental that y ^ 0. 
We can apply the above reasoning twice 



g(x,p)e h = 



g(x,p- h) -g(x,p) ^iiv 
— e h 



p&A, < 



e 2ir i y _ I 



E 
E 



g(x,p-h) -g(x,p) e 



2-Ki(p J r h)y 2tt i p y 

h — e s 



e 2?r j j/ _ J 



e 27ri 3/ _ I 



g(x,p- 2h) - 2g(x,p- h) + g(x,p) ^ lPy 

e 2w i y — I 



EO(h 2 ) 27vipy 
— - — e h 
g27r« y ^ 

p6A, < 

Let K C M be a bounded set such that 2irp) = if p ^ if, we 
can write 



TO 




.1 %)) C(z,y) e l ^" + « J - 



E2 TV % p y 
f{x,27rp)e~^~dxdy, 

pGhZHK, < 

and using the above calculation for g = f (note that derivatives of / 
with respect to p are well controlled) , we have that 



\I h 2 (e)\<0(tf) £ / / 



1 - /?%) 



e 2ni V — I 



C(x,y) ei s ^ 



dxdy 



0(h) / / dxdy dp 
Jk Jt J 



since 



1 - rf{y) 



D 2iri y J 



C(x, y)e^ s ^ y) 



< C, for some constant C. 



□ 

We now need to estimate 1^, for which we need to consider two cases 
la) and lb) depending on whether the support of / intersects or not 
the energy level H(P). 



Case la) Suppose that the support of / does intersect the energy 
level H(P). 

For each p we define a function 
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J j J j 

hence /f(e) = -^(p), where K is the p-projection of the sup- 

port of /. 

By the hypothesis of case I we are considering only points p^ in < 
P < Pmax- I n order to use the stationary phase method to estimate 
I zip), we need to analyze the points (x,y) such that VSp(x,y) = 0, 
it was shown that there are 4 points satisfying this, they are (xi(p), 0) 
and (0,2xj(p)), for % — 1,2. We need to calculate D 2 S(x,y) for these 
4 points. For (xj(p),0) we have that 

«^(<(4)) << f )) )' 

as 27rp 7^ P + (0), we have that Xi{p) 7^ 0. Hence ud(xj(p)) 7^ 0. Also 
V|det(^(^(p),0))| = |<(x t (p))| 7^ 0. 
For (0, 2xi(p)) we obtain 



D 2 S(0,2xiip)) 



v»(2xi(p)) 

<(2xi(p)) ; ' 



i = 1,2. 

To have Vo(2a;i(p)) ^ we need that 2xi{p) 7^ 0, this happens be- 
cause 2717) 7^ p + (0), and that = 2xj(j3) 7^ ±1, but these points are 
avoided because 7] € (1) = r] e (— 1) = 0. Hence y^jdet (D 2 S(0, 2xi(p)))\ = 
|<(2^(p))|^0. 

Now we use the stationary phase method, and note that Sp(xi(p), 0) = 
0, to obtain that 

= E uT^yTl /0*(p), 2^) C(^(P), 0) [1 + 0(v^)] + 

+ E rwW^^)) /(°> C(0, 2x. i (p))et^(°' 2 ^tf))[l+0(v^)]. 
Therefore 

#(e) = Jf + j£(e), 
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where 

2nh 



%K(x,(p))| 



and 



P&(hZ+£)nK,< i= 



4(e) 



E E 



2tt/i 



^ K(2xi(p))| 



• "rf (2x,GS)) /(0, 27rp) C(0, 2x l {p))e^ 2x ^\l + 0(\4)]" 

Remember that p + {xi(p)) = 27rj3, % — 1,2 (see Remark just after 
( I17jl ). For each pelwe define Xi(p), i — 1, 2 be the points satisfying 



(19) P + (^(p)) = 27rp. 

Therefore, when /i->0we have the following estimate 

f(xi(p),2np) t( Xi (p),0) 

and 



Jf(e)~27r^ / 



dp 



4(e)~27rJ2 



i=l 



rj e (2xi(p))f(0,27rp) C(0,2xi(p)) e^°' 2xi ^ 



K 



K'(2*i(p))| 



dp. 



The trouble here is that Vq(2xj(p)) is small if Xi(p) is near to 0. Let 
us estimate VQ(2xi(p)). 

We have Uo(2xj(p)) = ^o'(O) 2xj(p) + 0(xi(p) 2 ), we know that 



p + (x t (p)) = y/2[H(P) - V(x l (p))} and p max = y/2[H(P) - V 
Therefore 

p 2 max - (27vp) 2 = V( Xi (p)) - V(0) = Xi (p) V(0) + 0(^(p) 3 ), 
this implies 

Cxi(p) 2 = 2p max (p max - 2np) - (p max - 2irp) 2 + 0(x;(p) 3 ), 

then 

Xi(p) 2 (C + 0(xi(p)) = {p max - 2np)(D - (p max - 27rp)), 

hence 

*i(p) = (*w - 2ixp) [ C + 0(Xt(p)) 
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For e > small enough, if < e we have that \p ma x — 2np\ < e 2 C. 

And 



<)'l(2xi(p))\ « Cy/pmax ~ 2np + 0(p max - 2np) 



y/Pmax - 2np (C + 0(^p max - 2irp ) « C^Jp max - 2-np . 



Now we estimate J% (e) 



M 



pehznK,< i=i IWWJJI 



^(2xi(p)) 



2-E 



i=l 



M 

7j e (2z i (p))-=_= 

if yVmax ~ 



2 /• 

dp < M 2tt V / 

i=i M*: 



1 



y/Pmax ~ 2l\p 



dp, 



where C(e) = {p | < p ma x — 2ixp < Ce 2 }. The contribution over the 
variable p is like dp, which gives a factor of a/ Ce 2 . We point 

out that in the Riemann sum above (which does not contain the con- 
tribution of pd) the error of approximation for the integral J 6 dp 
is of order h. 

Finally, for each h fixed we consider p^ = Pd(h) a s before, and we 
analyze the term 

f(x,27rp h ) ((x,y) e^ s ^ x ^ dxdy. 

T JT 

Let 5 > be fixed, we define the set Bs = {(x,y);\x\ < 5, y G 
[— |, |]} U {(x,y); 1 2/ ] < 5, x e [— |, |]}. Then, as the stationary points 
are (xi(p),0) and (0,2xj(p)), for i — 1,2, we see that this points are 
in Hence | VSp h (x, y)\ > K(5) in J3|, where K(5) is uniformly 
bounded by below as h — > 0, for fixed 5. Therefore we have 





T JT 



< 



f(x,2irp h ) C(x,y) eh s ^ x >^ dxdy 
f(x,2irp h ) ((x,y) e^ Sf hi x >y) dxdy 



< 



+ 



+ 



f(x,2Trp h ) ((x,y) eh^h^^ dxdy 



Note that 



f{x,2-Kp h ) ((x,y) eh s Ph( x ^ dxdy 



< C5 
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And 



f(x,2irp h ) C(x,y) eh s vh( x >y) dxdy 



h f e iSf(x, y ) f&, ^V) C(x, y) n(x, y) VSpjx, y) ^ 
i .latnc\ \VSp(x, y)\ 2 



am 
h 



e iSp(x,y) y 



1 J B 



f(x, 2np) ((x, y) VS^x, y) 
\VS p (x,yW 



dxdy 



Hence 



f(x,2irp h ) ((x,y) e^ {x > y) dxdy 



< hK(6). 



Thus sending h — > 0, and then 5 — > we see that this term vanishes. 



Now we are able to estimate If(h). 

Remember that ((x,y) = e [-Mx+%)-"i( x -%)+ 2 Mx)} _ \y e will show in 
the last section that v\{x) = |ln(j> + (x)), hence ((xi(p),0) = p +^ X ( p )) 
and by the definition of v , v'^Xiip)) = (p + )'(xi(p)). Therefore 



J ! ^e)^2nJ2 J 



■fj \(p + )'(Xi(p))\p+(Xi(p)) 



dp. 



Above we are approximating an integral by a Riemann sum. 
Hence, as e, 5 > are arbitrary we get that 

A ] h' \(p + )'^p))\p + Mp)) P ' 



It is well known, via action angles variables (see |PR] |Ar] ) . that the 
projected Mather measure is given by the density 

dH(P) 



bp(x) = b(x) 



d 2 G, 



dP 



dxdP Pp(x) J 

were Gp(x) = j\ pp(s)ds. 

Using the change of coordinates Xi in each branch we get that for 
any given function g(x) we have 

f \ht \j f 9( x i{p))b{ x i(p)) , 
g{x) b(x)dx = / — f ~ ^ If — ^—^ — dp 



{p + )'(xi(p)) 
This describes the limit measure in this case. 
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We postpone the analysis of case lb) as it will use certain results 
which will be established in case II). 

Case II) Suppose that the support of / is contained in T x [c, d], with 
c,d such that [c,d] D [PmimPmaxl = 0- Hence, by item (ii) above, for 
each p £ KL + ^ such that l-np £ [c, d] we have V5*p(x, y) 7^ for all 
(x,y) £ T x T. Also, we note that |VS*p(x, y)\ 2 is uniformly bounded 
away from in all variables. 

We write 

i f (h)= Yl Mp), 

where 

hip) - I I f(x,27rp)((x,y)ei s ^ dxdy. 
Jt Jt 

Using the stationary phase method when TxT contains no stationary 
points (see |Mi] sec. 5.7, or Theorem 7.7.1 |Hoj ). we get 

I h (p) = - [ c is*(x,y) f( x > 27l P) C(x, y) n(x, y) T VSp(x, y) rfg 



h 



i Jd(TxT) \VSp(x,y)\ 
f(x, 2-np) C(x, y) VS p (x, y) 



e iS p (x,y) y 

1 ./TxT 



dxdy =: F + G, 



|V^(x,y)| 2 

where n is an exterior normal unit vector and dS denotes an element 
of surface area on the boundary d(T x T). 

Note that the integrand is not 1-periodic in y so the boundary term 
does not vanish. 

To analyze F, note that 



h 



f = -i e i s ^ x >y) f( x > y ^ n ( x > y) vs p( x ' y) dS 

i Jd(jxT) \VSp(x,y)\ 2 

h H i 3 . (x uf(^p)C(x,-i) d -§(x,-^) 

h p( ' 2> ; : tttt; dx+ 



i 7-1 ~ \VS p (x,-^ 2 



2 ■ 



\VS p (iy)\ 2 dy+ 



2 



h n ± s -(x i)f( x ' 27r p) c{x,\)^{x,\) 

J / ph J P\ J '>2' 2 d,T — 
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As the functions Sp, f and £ are periodic in the variable x, the sum 
of second and the fourth term on the right hand side of the equality 
above is zero. Also, by the periodicity in x of the functions involved, 
the first and the third terms are equal to 0. 

Hence F = 0. 

Let us now analyze G, define 

, , V7 \f(x,2irp)C(x,y)VSp(x,y)' 

gp{x,y) := V — — 

[ \VSp(x,y)\ 2 

we can repeat the calculation and get 



G 



ASf(x,y) 9p(x, y) n(x, y) T VSp(x, y) 
\VSp(x,y)\* 



d(TxT) 



TxT 



gp(x,y) VSp(x,y) 



\VSp(x,y)\- 



dxdy. 



Then, as If(h) = F + G, we have 



p£hZ 



If(h) = h 2 



eh 



s 6 (x,y) 9p(x, y) n(x, y) T VSp(x, y) ^ 



a(TxT) 



\VS,(x,y)\- 



E 



TxT 



gp(x,y) VSp(x,y) 
\VSp(x,y)\ 2 



dxdy =: G\ + G 2 - 



We will show that \Gi\ = 0(h), by similar arguments we can show 
that \G 2 \ = 0(h). As \ e i s ^\ = 1, we have 

gp(x,y) n(x,y) T V Sp(x,y) 



\GA < h 



h 



* E 



pehi 



9(TxT) 



\VSp(x,y)\< 




g p (x, y) n(x, y) T VS p (x, y) 



dSdp 



<9(TxT) 



dS 



0(h), 



\VS p (x,y)\ 2 

the last equality is true because p G [c, d] , all functions involved are 
sufficiently smooth and \V S p (x,y)\ 2 is uniformly bounded away from 



0. 



Case lb) Suppose that the support of / does not intersect the energy 
level H(P). 
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We have already shown that /2(e) = 0(h), hence we just need to 
analyze 



^00 = E (v e (y))f(x,^P)t(x,y)e mx ' v) dxdy. 



We will suppose that supp(/) C [a, b] x [c, d] = A, where A is such 
that A n | x G [-§, |]} = 0. 

For each p fixed the critical points (xi(p), 0) for S^x, y) do not con- 
tribute in the integral, because f(xi(p), 2irp) = 0. Then or we have the 
contribution of the points (0, 2xi(p)), that we have already shown that 
is of order O(e), or in A there are not stationary points. In this last 
case, there exists Kp > such that \V Sp(x,y)\ 2 > Kp if f(x, 2np) ^ 0. 
As p lies in a compact set there exists K > 0, uniformly in p, such that 
\VSp(x,y)\ 2 >Kiif(x,2np)^0. 

Also, note that the function r] e (y) is smooth and periodic in y. There- 
fore we can apply the same argument used in case II) to prove that 

\iH*)\ = o{h). 

Proof, (of Theorem [TJ By the analysis above we see that the limit 
measure is supported in the graph {(x, p + (x)) \ x G [— \, |]}. 

Item a) follows by case la), item b) follows by case lb) and case II) 
above. □ 



3. The asymptotic expansion of v h 



Let Vh and uj* be T-periodic solutions of 

-^ + \\P + v' h \ 2 + V = H h (P), 

and, 

h ^f + l\P + v* h '\ 2 + V = H h (P). 

We will provide rigorous asymptotic expansions for Vh and v* h , as 
required in the previous sections. We will only present the proofs for 
Vh as the analysis of is completely similar. 

For a given N we will construct functions Vi, numbers Hi, i = 
0, N, and formal expansions 

i)h = v + h vi + h 2 v 2 + ... + h N v N , 
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and 

= H + hHx + h 2 H 2 + ... + h N H N . 

We will show that 

H h = H h {P) = H^ + 0{h N+1 ), 

and 

v h (x) - v h (x h ) = (x) - v% (x h ) + 0(h^). 

In order to get these estimates we will derive differential equations 
which define and Hi in a suitable way. In fact, vq and H are solutions 
of the limit problem 

H = H(v' (x),x) = V(x) + \(P + v' (x)) 2 . * (1) 

The remaining equations are obtained by Taylor expansions in formal 
powers of h. For instance, the following equations will be 

#1= -^ + HM,*W= -Y + (P + v> )v' 1: *(2) 

and 

v" (v 1 ) 2 
Hi= -f + HM,x)4 + ±f-. *(3) 

All approximations in this section are 0(h~^~) in the uniform topol- 
ogy- 

From *(1) we have that v = (p and H = H(P), where <f){x) = 
f-i Pp(s)ds — P (x+ |). Using this last expression in *(2), we have an 

equation for v 1 and H\. The value Hi is obtained by the compatibility 
condition: 

Hi = - J Vq(1<to, 

where o"o is the projected Mather probability for h — 0, which satisfies 
div {(P + v' )a ) = 0. 

The property div ((P + v ) a"o) = follows from the fact that the 
Mather measure \i has support on the graph (x,P + v (x)), and also 
satisfies the holonomic condition J vD x tp(x)d^ = 0. This implies for 
the projected measure o"o that J(P + v (x))(p x dao = (see |Evlj . |BG] 
or [EOT] ). 

Then, once determined Hi, the function vi is uniquely determined 
(up to additive constant) because (P + v' ) is never zero. 
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Remember that v'q = (pp)' and that the density of the projected 

d~H{P) 

Mather measure is given by bpix) = ? f : . Hence 

— dH(P) ft (pt(x))' , d~H(P) , , ,,i 



and the equation *(2) becomes 

= -M^ 



Therefore 

= -ln(p£(x)). 

Using the function i>i in *(3) we have an equation for v 2 and i?2- 
There exists a unique number H 2 , given by the compatibility condition 

J \2 



to which corresponds a unique (up to constant) v 2 . Again, the solv- 
ability is ensured by the condition (P + v' ) ^ 0. 

In this way we obtain, inductively, all the i>-, and all the Hj, j = 
0, 1, N. Note that the functions Vj are smooth, because V is smooth 
and consequently p + also is. 

Proposition 2. For a// N and for all x 

(20) - ^im + 1 | P + + = + o^+i). 

Remark. As the proof below shows, because the expansion vj^ is a 
finite sum of smooth functions, the error term in the right hand-side of 
$20]) is 0(h N+1 ) in the C k topology for all k > 0. 
Proof. We will show the claim for N = 2, as all other cases are 
analogous. We want to show that 



h(v 



2\ni 



+ H((ti)'(x),x)-HZ = 0(h 



2 

In order to do so, we expand, for each fixed x, H((v%)'(x),x) in 
Taylor series: 

H((v 2 h )'(x),x) = 
h 2 

= H{v' ,x) + hHpiv^x)^ + hv' 2 ) + —Hppivix)^ + hv' 2 ) 2 + 0(h 3 ). 
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Hence, by definition of Vi,Hi, i = 0, 1, 2 we have 

_mm +H{mx)tX) _ 6 2 = 

= -7?« + hv'i + h 2 v' 2 r ) + H(v' , x) + hH p (v' , x)(v[ + hv' 2 )+ 

+yi?RpK, x)(v[ + hv' 2 ) 2 + 0(h 3 ) - (H + ^ + h 2 H 2 ) = 
= h 3 [-v' 2 ' + H pp (v' ,x)(2v[v' 2 + h(v' 2 ) 2 )} + 0(h 3 ) = 0(h 3 ). 

□ 

Proposition 3. \H h - Hjf\ = 0(h N+1 ) 
Proof. 

From |Gomlj . section 3, we have that 

H h (P) < sup { - Whp°l + H((v»)>(x),x) }<H? + O(h^). 

Therefore, H h - Hf? < 0(h N+1 ). 

As H(p, x) is convex on p, we have that 

H{v' h {x),x) - H((v»)>(x),x) - \ «(*) - «)"(*)) + 0(h N+1 ) > 

H p (($»)'(x),x) «(*) - (v»)'(x) ) - \ K(x) - (t>£7(*)) + 0(h N+1 ). 
Consider a point xq where (vh(x) — (x)) has a maximum. Then, 

K(*o) - (W(^o))_= and h K(x ) - (Of )"(*o)) < 0. 
From the above H h — H™ > 0{h N+1 ). □ 

Lemma 2. There exists a periodic non negative probability density cxf 
that is a solution of the equation 



h 



where g^(x) = H p ((v^)'(x),x). Also, there exist constants k, K > 



such that k < cj^(x) < K . 



Proof. Note that g% — P + v' + hv[ + h 2 v 2 is positive and bounded 
uniformly in h. For convenience we will drop the N. By integration, 
we get that the equality for all x 

{gh{x)a h (x))' + hah = 0, 
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is equivalent to 

g h (x) a h (x) + ha' h {x) = c h , 

for some suitable constant Ch- 

We multiply the above expression in each side by the integrating 
factor oth(x) = e~h~ , where G(x) = J* gh{s)ds. 

In this way, 

h {a h a h )' = ha' h (x)a h (x) + ha h (x)a' h (x) = 
&h{x) g h {x) cr h (x) + ha h (x)a' h (x) = c h a h (x). 

Therefore, 



ha h (x)a h (x) = c h a h (s) ds + C h , 
Jo 

or, equivalently, 

°h(x) = ^ . 

The constants Ch and Ch are adjusted so that ah is a periodic prob- 
ability density. 

We claim that the constants Ch and Ch are bounded. Note first that 

g h {x) a h {x) dx = / [g(x) a h (x) + ha' h (x)] dx = c h . 
Jo 

As gh is positive and bounded uniformly on h, and ah is a density, 
we get that Ch is bounded. In fact, Ch converges to a constant, that we 
will denote by c. 

By Laplace method, see [01] (chapter 3, section 7) |Mij . as G is 
monotonous increasing, we have 

Jo e h ds 1 



Furthermore, the limit is uniformly bounded on h. Now we will show 
that Ch is uniformly bounded in h. 
As 

a h (0) = a h {l), 

we get 

Ch C h C h + c h ds 
n e >>■ n e h n 
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Consequently, 

„ „ am c h Jeh ds 
C h = C h e * + J0 M h. 

e >>■ h 

From this 

am , Ch L e h as 
e ^ /z 

In this way, — > 0, as /i — > 0. Therefore for any x 

Ch C cc/i(s) (is + C/j c/j /!f e - ^" ds 
hmWx) = Jim ^ ^ = lim u — 

> 



P + v' (x) 

This shows that ah is bounded and bounded away from zero. □ 
Proposition 4. 

(21) Wh-{v N h )'\ L , = 0{h^). 

Proof. 

Note that from the last proposition 

0(h N+1 ) = 

H((v»)>(x),x) - k - H(v' h (x), x)) + ^ 

From the strict convexity of the Hamiltonian, and because H(p, x) 
is quadratic in p we have the following inequality 

H(v> h (x),x) - If («)'(*),*) - \ «'(*) - «)"(*)) > 



ajf (x)dx. 



H p ((v»)\x),x)(v' h (x)-(v»)'(x)) + l(v' h (x)-(v»)>(x)) 2 



xw-K v rw). 

Let a 1 ^ be as in the previous Lemma. We claim that for any periodic 
function w we have 

(22) J (H p ((v»)'(x),x)w'(x) - ^Vf(x) = 0. 

This follows from integration by part and from the last lemma, with 
9h( x ) = Hp{(Vh)'{ x )i x )i f° r a ^ x we have 

[ ( H p ((e»y(x),x) af^x) y - &otf )"(*) 1 = o. 
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It follows from ( 1221) applied to w = Vh — that 

/ «(*) - {v»)'{x)Ya»{x)dx = 0(h N+1 ). 
As is bounded and bounded away from zero we get 
/ \v , h {x)-ffi)\x)\ 2 dx = 0{h N+1 ), 

and, so 

\\v' h -m\\ L ^o(h^). 

□ 

Proposition 5. 

(23) \v h (x) - v h (y) - (v»(x) - v»(y))\ = O(h^). 

Remark. The expansion in (123|) holds only in the uniform topology, 
as Vh is in general only continuous. Nevertheless the function (v^ is 
smooth. 

Proof. Define gh(x) = Vh(x) — v^(x), this function is T-periodic. 
Using Morrey's Inequality with n = l,p — 2, we have that 

\9h(x) -g h (y)\ < C\\g' h \\ L2m , 

for all x, y G T. By proposition H] we get the result. 

□ 

The analysis of v\ is similar to the case we just described. 

We would like to thanks the referee for many valuable comments and 
for a very careful reading of the manuscript. 
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